First-principles study of the spin-lattice coulpling in spin frustrated DyMn2 05 
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The lattice dynamic properties and spin-phonon coupling in DyMn2 05 are studied by using the 
density-functional theory. The calculated phonon frequencies are in very good agreement with 
experiments. We then compare the phonon modes calculated from different spin configurations. 
The results show that the phonon frequencies change substantially in different spin configurations, 
suggesting that the spin-phonon coupling in this material is very strong. Especially, the short 
range spin ordering has drastic effects on the highest Raman and IR phonon modes that might be 
responsible for the observed phonon anomalies near and above the magnetic phase transitions. 

PACS numbers: 75.25.+Z, 77.80.-e, 63.20.-e 



I. INTRODUCTION 

Manganese oxides including RMnOsjii^ and RMn2 05 
(R=Tb, Dy, Ho, etc.}^i^i^ are a very special class of im- 
proper ferroelectrics that display strong magnetoelectric 
(ME) effects, such as the "colossal magnetodielectric" 
(CMD) effects and magneto-polarization-flop efFectsj^'^'^ 
These materials have attracted great attention,— i^i^iiiS be- 
cause the strong ME effects can be utilized to tune the 
electric properties of the materials via applied magnetic 
filed, or vice versa, and therefore have great potential for 
future multifunctional device applications. The macro- 
scopic ME effects have their microscopic origin from the 
interplay between the lattice degree of freedom and spin 
degree freedom. Neutron scatteringjii as well as first- 
principles calculationsi^ii^ shows that the improper fer- 
roelectricity in TbMn2 05 is driven by the non-central 
symmetric magnetic ordering, which, by coupling to the 
lattice, lower the crystal symmetry] -'^^'^^ The spin-lattice 
coupling can also modify the lattice dynamitic prop- 
erties. Indeed, recent temperature dependent Raman 
measurements^'^ show anomalous phonon shift at T* ~ 
1.5 Tjv and near the magnetic phase transition tempera- 
ture T^r in DyMn2 05. Similar anomalies are also ob- 
served for the IR modesJ^ Further B-field dependent 
measurements show that the IR modes are very sensi- 
tive to local magnetic structure^^ which suggest that the 
spin-phonon coupling in this compound is very strong. 
However, the experiments give only the overall effects, 
and the detailed mechanism of the spin-phonon coupling 
is still lack. 

First principles method has been applied successfully 
in studying spin-phonon coupling in geometrically frus- 
trated spinel ZnCr204J^ It is a challenge problem to 
study the spin-phonon coupling in the spin frustrated sys- 
tems, such as RMn2 05, because in these materials, the 
magnetic interactions are of different magnitudes (e.g., 
\Ji\i \Jb\ ^ l-^sDr^*^ and the spin frustration leads to 
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complicated spin-spin correlations as functions of tem- 
perature and magnetic field. Therefore, the spin will 
develop different short range ordering at different tem- 
perature above the magnetic phase transition, and make 
it hard to identify the origin of the phonon anomalies.— 
In this work, we study the spin-phonon coupling in 
DyMn205, to shed some light on the observed phonon 
anomalies. We show that the short range spin correla- 
tion and local spin structure have significant effects to 
the phonons frequencies in this material. 



II. METHODOLOGY 

The first-principles density-functional (DF) calcula- 
tions are done by using the Vienna ah initio Simula- 
tions Package (VASP).— — We use a spin-polarized gen- 
eralized gradient approximation (GGA), with Perdew- 
Burke-Ernzerhof functional. ^° A plane-wave basis and 
projector augmented-wave (PAW) pseudopotentialsSi are 
used, with Mn 3p3(i4s, and Dy 5p5(i6s electrons treated 
self-consistently. A 500 eV plane-wave cutoff results in 
good convergence of the total energies. We relax the 
structure until the changes of total energy in the self- 
consistent calculations are less than 10~^ eV, and the 
remaining forces are less than 1 mcV/A. DyMn205 un- 
dergo several magnetic phase transitions^ at low tem- 
perature. We are interested in the phonon properties in 
the paramagnetic (PM) and anti-ferromagnetic (AFM) 
phase. To accommodate the magnetic structure, we use 
a 2 X 1 X 1 supercell.^^ For the supercell we used, a 1 x 2 x 4 
Monkhorst-Pack k-points mesh converges very well the 
results. 



III. RESULTS AND DISCUSSION 

The fully relaxed DyMn2 05 is an AFM insulator, and 
has same lattice structure to that of TbMn2 05,^'?i?-- both 
are of P62im symmetry, as Tb and Dy are isovalent. The 
lattice structure is distorted from a Pham high [H) sym- 
metry structure, due to the spin-lattice coupling. ^^ The 



TABLE I: Calculated phonon frequencies of ground state structure L. The exp eriment vaules for Raman phonons are extracted 
from Ref. [l4l . whereas the IR phonons frequencies are extracted from Ref. [la . 
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calculated lattice constants are a=7.270 A, 5=8.518 A 
and c=5.600 A, respectively, and are in very good agree- 
ment in experimental values^ (7.285, 8.487 and 5.668 
A respectively). The lattice constants are somewhat 
smaller than those of TbMnzOs J^i22 Like ThMzOg^^i^ 
DyMn2 05 has two energetically degenerate lattice (and 
magnetic) structures, L and i?, in which Mn*+ form an 
AFM square lattice in the ah plane, whereas Mn'^+ cou- 
ples to Mn^"*" either antiferromagnetically via J4 along a 
axis or with alternating sign via J3 along b axis. Mn'^+ 
ions in two connected pyramids also couple antiferromag- 
netically through J5. Here, we adopt the notations J3, 
J4 and J5 from Ref. Q, and define the J3 to be the Mn^+- 
Mn'^"'" superexchange interaction through pyramidal base 
corners, and J4 the superexchange interaction through 
the pyramidal apex. 

We first calculate all zone-center optical phonon fre- 
quencies of the fully relaxed structure (L or R) via a 
frozen-phonon techniquCfi^ and the results are listed in 
Table HI Symmetry analysisi^ shows that for the high 
symmetry (Pbam) structure H, phonons belong to 8 ir- 
reducible representations (irreps), among them i?i„, i?2« 
and B^u are IR active, where Big, B2g, B^g and Ag are 
Raman active and A„ is silent. However, in the low sym- 
metry {Pb21m) structure {L or R), the Raman and IR 
modes couple to each other, and regroup to 4 irreps, i.e.. 



A-l — B2u 

Bi — Biu 



Ag; A2 = B 



2ff 



Au] B2 — B^u ffi B 



Iff' 



B^g. The couplings between irreps. are 



very small, therefore the results are given by their ma- 
jor symmetry characters. As we see from Table U the 
results are in excellent agreement with experiments i^^i^^ 
The calculated phonon frequencies are very close to those 
of TbMn205, because Tb and Dy are isovalent and the 
two materials have very similar lattice structures. The in- 
ner products between the corresponding phonon modest 
of the two compounds are close to unity, suggesting that 
they have very similar mode patterns. 

Experimentally, there are several anomalies found for 
both Raman modes and IR modes, near the magnetic 
phase transitions. For example, anomalous phonon shift 
for the Raman modes had been observed at T* ^^ 60-65 
K, and near the Neel temperature T^r J^ The highest Ag 
mode show a steep hardening upon cooling between T* 
and TAffii For the IR modes, it has been found near 60K, 
the infrared active modes soften and several modes soften 
substantially.^^ The phonon modes are also found to be 
very sensitive to the applied magnetic field. ^^ The ob- 
served phonon anomalous phonon shift strongly suggest 
that the spin-phonon coupling in this material is very 
strong. 

Theoretically, the spin-lattice coupling in this com- 
pound can be described via an effect Hamiltonianji^ 



J 



E{{ux}) = (^0 
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where u\ is the A-th zone-center phonon, and S^ is the 
magnetic moment of the i-th atom. Here, we consider 



only the magnetic moments of the Mn^+ and Mn^+ ions. 
Jij is the exchange interaction between two adjacent Mn 







FIG. 1: (Color online) The spin configurations used in calcu- 
lating phonon frequencies of different spin states. 



ions. Eq is the total energy of high symmetry structure 
H without spin-spin interactions, where the ground state 
spin configuration is determined by the first term, lox 
is the frequencies of the A-th phonon, calculated in the 
high symmetry PM state. The spin-lattice coupling have 
two effects: (i) The second term, which is linear in u, 
breaks the inversion symmetry of the system, and drive 
the structure to a low symmetry polarized state leading 
to coexistence of AFM and ferroelectricity. The coupling 
strength is proportional to dJij / du\. However, this term 
would not change the phonon frequencies (at least to the 
first order approximation), (ii) The phonon modes ujx 
(including frequencies and eigenvectors) are modified in 
the presence of spin-phonon interaction, due to the third 
term of Eq. ([T]) , where the spin-phonon coupling strength 
is d'^Jij/duxdup. Different spin configuration (SC) {8^} 
would therefore have different phonon frequencies. 

To investigate how the SCs change the phonons in 
DyMn205, we calculate the phonons frequencies of dif- 
ferent SCs including the high temperature PM state, the 
AFM state, and the ferromagnetic (FM) state. The spin- 
lattice coupling strength J" can be extracted by com- 
paring the force-constant matrices calculated from these 
different spin states.— [see Eq. ([T])]. In this work, we fo- 
cus on the B2„ and Ag modes, whereas other modes can 
be studied in similar way. To simplify the discussion, 



all following calculations are done in the high symmetry 
structure iJ, constructed via symmetrizing structures L 
and R according to the Pbam symmetry.— In reality, the 
lattice constants would be somewhat different at differ- 
ent SCs. In the present calculations, this effect is ignored 
and the lattice constants are fixed in the calculations. 

The phonon frequencies of PM, AFM (SC G in Fig. [1]) 
and FM (SC H in Fig. [T|) states are compared in Table HTl 
for the Ag modes and in Table IIIII for the 62^ modes. 
One may attempt to calculate the phonon frequencies of 
the PM state using spin unpolarized GGA, which is listed 
in the column under u-GGA. A quick look reveals that 
the phonon frequencies in this column are significantly 
lower than those of AFM and FM state. Especially, as 
shown in Table |TT1 the Ag irrep has a soft mode, caus- 
ing a Jahn- Teller distortion.— The phonon frequencies 
calculated from u-GGA is not a good approximation for 
the PM state, because the Mn ions have local magnetic 
moments even in the PM phase, although their direc- 
tions are distributed randomly. Alternatively, we calcu- 
late {u!)dis, which are the averaged phonon frequencies 
of several fully disordered SCs {A, B, C in Fig. [T] ) so 
all exchange interactions are canceled out. The averaged 
phonon frequencies {uj)dis are close to those of the AFM 
and FM states without soft phonons. 

However, the spins arc fully disordered only at very 
high temperature. At lower temperature, especially close 
to the magnetic phase transition, the spins are somehow 
correlated, and develop short range ordering. Since J4 
is the largest among all the exchange interactions,^- to 
compare with experiments, we consider the ordering of 
J4 interactions. We calculate {uj),j^, which are the aver- 
aged phonon frequencies of three SCs (D, E, F in Fig. [T]) 
in which J4 are fixed in the AFM configuration, whereas 
all other exchange interactions are canceled out. As we 
see from Table |lT] and Table IIIII the frequencies of most 
modes in different SCs differ by about 3-5 cm~^, which 
are of similar magnitude to the phonon shifts in the mag- 
netic field experiment. ^^ However, the high frequency 62^ 
(671 cm^^) and Ag (679 cm~^) modes change dramat- 
ically (~ 20 cm^^) for different SCs. Especially, {uj)dis 
and {uj)j^ are different by about 10 and 13 cm~^ for the 
Ag and B2U modes, respectively, though both belong to 
the PM phase. Interestingly, the two modes frequencies 
of {uj)j^ are close to the those in AFM state, and {uj)dis 
have the mean values of those of AFM and FM states. It 
is usually believed that GGA overestimates the exchange 
interactions, whereas including the on-site Coulomb cor- 
relation may improve this problem. To see how includ- 
ing the on-site Coulomb correlation will change the spin- 
phonon coupling in DyMn205, we have carried out GGA 
+ U calculations of the phonons in different SCs. The 
on-site Coulomb U of Mn ions has been taken to be a rea- 
sonable value 4.0 eV, whereas the exchange parameter j— 
0.8 eV is used. The results are also listed in Table HIl and 
Table Hill for the Ag and 62^ modes respectively. As we 
see, GGA + U significantly overestimates the frequencies 
of the high frequency modes when compared to experi- 



TABLE II: 
methods. 



The Ag modes calculated in high symmetry structure H of different spin configurations using GGA and GGA + U 
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{u)ji 
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{u):ji 


AFM 
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571.1 


667.7 


677.6 


679.0 


662.8 


746.1 


759.9 


759.0 


723.2 


563.2 


621.5 


625.8 


618.6 


629.7 


658.8 


661.3 


659.4 


653.4 


530.7 


617.7 


618.4 


613.5 


617.3 


634.2 


635.3 


631.0 


638.3 


463.3 


537.2 


540.9 


531.1 


547.0 


559.3 


561.1 


556.6 


562.2 


432.9 


490.6 


492.2 


491.3 


494.1 


505.4 


506.3 


506.0 


506.4 


409.2 


441.0 


442.8 


442.6 


441.4 


451.4 


452.7 


452.6 


450.1 


369.6 


414.3 


407.9 


407.0 


429.8 


438.2 


436.6 


436.0 


440.3 


249.6 


340.0 


337.6 


337.6 


340.3 


342.6 


342.3 


341.6 
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321.9 


329.2 
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331.7 
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331.9 
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334.8 
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236.2 


243.9 


248.3 


248.0 


247.5 


249.2 


102.9 


222.3 


220.2 


219.6 


225.9 


217.7 


216.9 


216.6 


220.1 


66.1 


136.8 


137.5 


136.0 


137.6 


141.1 


141.6 


140.9 


141.0 
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TABLE III: The B2u modes calculated in high symmetry structure H of different spin configurations using GGA and GGA -f 
U methods. 
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FM 
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586.4 


660.2 


673.3 


671.7 


649.6 


723.7 


736.6 


735.9 


705.8 


512.8 


627.2 


625.5 


626.1 


623.8 


659.1 


659.4 


659.1 


657.9 


507.2 


551.6 


552.2 


553.4 


552.3 


578.5 


578.7 


579.4 


577.9 


488.7 


534.0 


530.9 


535.3 


534.1 


566.6 


567.0 


568.3 


564.6 


431.2 


474.5 


473.5 


475.7 


476.4 


495.5 


494.5 


495.6 


498.1 


391.0 


439.5 


437.2 


439.2 


438.1 


461.8 


462.0 


462.2 


459.4 


381.3 


409.7 


410.0 


411.7 


406.0 


419.9 


420.2 


420.9 


418.0 


310.0 


353.6 


350.2 


350.8 


357.6 


365.7 


364.9 


365.2 


366.1 


274.3 


315.0 


314.3 


315.2 


314.6 


316.2 


315.8 


316.3 


316.6 


257.7 


272.7 


266.0 


267.9 


280.1 


283.8 


282.5 


282.5 


285.5 


211.1 


222.2 


221.8 


223.4 


220.8 


211.7 


212.0 


212.5 


212.7 


142.6 


162.0 


162.2 


162.0 


161.4 


165.6 


165.7 


165.5 


165.7 


139.2 


156.8 


156.9 


157.0 


157.3 


163.3 


163.5 


163.7 


162.7 


93.9 


98.4 


98.6 


99.2 


96.8 


94.0 


94.3 


94.6 


93.4 



ments. However, the frequency difference between differ- 
ent SCs are almost remain the same to those without [/, 
suggesting that the spin-phonon couphng constants do 
not change much by inchiding the on-site Coulomb cor- 
relation. In the following discussion, we therefore focus 
only the GGA results. 

To understand the results, we analyze the vibrational 
pattern of the two modes, shown in Fig. [2l For the 
high-frequency Ag mode, about 90% of the vibration is 
associated with O3 (the oxygens that connect Mn'^+ and 
Mn'*+ along the a axis, and convey J4 interaction). The 
O3 atoms vibrate in the ab plane, with an angle of 21.4° 
to a axis. This mode also has small components of Mn'^+ 
vibrating in the a-axis, and Mn^+ motion in the c-axis. 
Therefore J5 may also affect the Ag mode (but consider- 
ably smaller than J4). In the high-frequency B2U mode, 
more than 70% of the contribution comes from O3, which 
also vibrate in the ah plane, with an angle of 19.9° to a 



axis. The rest contribution includes the motion of Mn^+ 
along a axis. If the spin exchange interaction is local, the 
motion the O3 atoms would only tune the J4 interactions. 
Therefore the J4 interaction plays an essential role of the 
spin-phonon coupling in these two modes. Now it is easy 
to understand the phonon shifts of the two modes in dif- 
ferent SCs. According to Eq. ([1]), the frequency shifts 
of the highest Ag and B2U modes are (S3 • S4) d^J^jdv?^, 
where S3, S4 are the spin vectors of Mn^+ and Mn'*+ 
ions associated with J4, respectively. For the fully disor- 
dered SC, (S3 • S4)=0. In the J4 short range orderd state 
and in AFM state (S3 • S4)=-S'3S'4, whereas in FM state, 
(S3 • S4)=S'3iS'4. iS'3=2.3 [iB and 54=1.64 /i^, are the lo- 
cal magnetic moments. By comparing the force-constant 
matrices of different SCs, we estimate d^Ji/du\ ~ 0.120 
meV/(A-^B)^ for the two high frequency modes. As dis- 
cussed in previous section, this value remain unchanged 
when including the on-site Coulomb interactions. 




oDy 



o Mn 



oO 



FIG. 2: (Color online) The mode pattern of the highest (a) 
Ag mode and (b) B2u mode. 



Experimentally, such large phonon hardening due to J4 
ordering was not directly observed during the tempera- 
ture dependent measuremen t^'^i^^ for two reasons. First, 
the Ji ordering develop gradually with decreasing of the 
temperature. A Monte Carlo simulations^ show that the 
average of (S3 • S4) increases from 80% to 90% of its 
maximum value when temperature lower from 1.5 Tjv to 
T^v. Second, the phonon frequencies hardening due to J4 
ordering is accompanied by the anharmonic effects j^ii^ 
and it is hard to isolate the anharmonic effects and spin- 
phonon effects in experiments. The lower frequencies 
modes involve collective motion of many atoms. There- 
fore, the phonon frequency shifts due to spin-phonon cou- 



pling are a consequence of the competition of many J"s, 
and are much smaller than the high frequency modes. 

Below Neel temperature, Tn, the material is in the 
long range ordered AFM state. We then compare the 
phonon frequencies of short range ordered state (w)j^ 
to those of AFM state. As we see, while the highest Ag 
mode hardens by 2 cm~^, and the highest B2„ mode softs 
by 2 cm~^, consistent with experiments^ I ~'^^ The fre- 
quencies difference between (w)j^ and AFM might come 
from the J5 interactions. 



IV. SUMMARY 

To summarize, we have investigated the spin-phonon 
coupling in a spin frustrated DyMn2 05 system via first- 
principles calculations. We compare the phonon modes 
calculated from different spin configurations. The results 
show that the short range spin ordering can dramatically 
change the phonon frequencies in this compound, and 
might be responsible for the observed phonon anomalies 
near and above the magnetic phase transitions. 
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